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^ \ Abstract 
> '■ 

■ The hadronic decay B J/t/jK* is analyzed within the framework of QCD factorization. The 

^ ' spin ampUtudes Aq, A^^ and Aj_ in the transversity basis and their relative phases are studied 

^ using various different form-factor models for B — K* transition. The effective parameters 02 for 

helicity h = 0,+,— states receive different nonfactorizable contributions and hence they are helicity 
dependent, contrary to naive factorization where 02 are universal and polarization independent. 
QCD factorization breaks down even at the twist-2 level for transverse hard spectator interactions. 
Although a nontrivial strong phase for the j4|| amplitude can be achieved by adjusting the phase of 
an infrared divergent contribution, the present QCD factorization calculation cannot say anything 
definite about the phase (/)||. Unlike B — > J/tpK decays, the longitudinal parameter 02 for B — > 
J/ipK* does not receive twist-3 corrections and is not large enough to account for the observed 
branching ratio and the fraction of longitudinal polarization. Possible enhancement mechanisms 
for 0° are discussed. 
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I. INTRODUCTION 



It has been well known that the factorization approach (naive or generalized) fails to 
explain the production ratio R = B{B J /ipK*)/B{B J/ipK) and the fraction of 
longitudinal polarization Tl/T m. B ^ J /ipK* decay. We consider two representative form- 
factor models for B — K{K*) transitions, the Ball-Braun (BB) model based on the light-cone 
sum rule (LCSR) analysis [0] and the Melikhov-Stech (MS) model based on the constituent 
quark picture. Both are consistent with the lattice calculation at large g^, the constraint from 
B (pK* at lower and the constraint from heavy quark symmetry on the dependence 
of heavy-light transition (see Sec. IV for more details). We see from Table I that in general 
the predicted longitudinal polarization is too small, whereas the production ratio is too large. 

TABLE I. The ratio of vector meson to pseudoscalar production R and the longitudinal polar- 
ization fraction Tl/T in B ^ J/ijjK^*'^ decays calculated in two representative form-factor models 
using the factorization hypothesis. 



BB 



MS 



Experiments 



CDF II 



CLEO § 



BaBar || 



Belle § 



R 



3.40 
0.47 



3.11 
0.46 



1.53 ±0.32 
0.61 ±0.14 



1.45 ±0.26 
0.52 ±0.08 



1.38 ±0.11 
0.60 ± 0.04 



1.43 ±0.13 
0.60 ± 0.05 



This is understandable because the parameter 02, which governs B — » J/iIjK{K*) decays, 
is assumed to be universal according to the factorization hypothesis, namely a2{J /ipK*) = 

a2{J /ipK) where /i = 0, ±, — refer to the helicity states 00, ±± and respectively. In the 

above-mentioned form-factor models, one has = 5.98, /i+ = 6.23 and /i_ = 043 (in units 
of GeV^) in the BB model and Hq = 5.47, /i+ = 5.92 and h_ = 0.73 in the MS model, where 
hi are the helicity amplitudes given by 



2mK 



rriB + rriK* 



'1.1] 



It is obvious that > Hq ^ h_. Therefore, under naive factorization Tl/T ^ 
(a2/io)^/[('^2^o)^ + (0-2^^+)^] = hl/{hl + h\) < 1/2 and R is expected to be greater than 
unity due to three polarization states for J/ipK*. These two problems will be circumvented 
if nonfactorized terms contribute differently to each helicity amplitude and to different decay 
modes so that al{J/ipK*) > atiJ/tpK*) ^ a^iJ/tpK*) and a2{J/ipK) > a'l^J / ip K*) . In 
other words, the present data imply that the effective parameter should be non-universal 
and polarization dependent. Recently two of us have analyzed charmless B — > VV decays 
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within the framework of QCD factorization 0. We show that, contrary to phenomenological 
generahzed factorization, nonfactorizable corrections to each partial-wave or hehcity amph- 
tude are not the same; the effective parameters Oj vary for different hehcity amphtudes. The 
purpose of the present paper is to study the nonfactorizable effects in B ^ J/ipK* decay 
within the same framework of QCD factorization. 

The decays B J/iIjK{K*) are of great interest as experimentally only a few color 
suppressed modes in hadronic B decays have been measured so far. The recent measurement 
by BaBar ^ has confirmed the earlier CDF observation that there is a nontrivial strong 
phase difference between polarized amplitudes, indicating final-state interactions. However, 
no such evidence is seen by CLEO |^ and more recently by Belle 0. It is interesting to 
check if the current approach for B hadronic decays predicts a departure from factorization. 
Therefore, the measurements of various helicity amplitudes in B ^ J/ipK* decays will 
provide a nice ground for testing factorization and differentiating various theory approaches 
in which the calculated nonfactorizable terms have real and imaginary parts. 

It is known that in the QCD factorization approach the coefficient 02 is severely sup- 
pressed in the absence of hard spectator interactions. It has been shown in p| that \a2\ in 
B — > J/ipK is of order 0.11 to the leading twist order, to be compared with the experimental 
value of order 0.25. The twist-3 effect in hard spectator interactions will enhance 02 to the 
value of 0.19to;i2- We shall see later that, contrary to the J/ipK case, in 5 ^ J /ipK* does 
not receive twist-3 contributions and it is dominated by twist-2 hard spectator interactions. 

The layout of the present paper is as follows. In Sec. II we first outline the necessary 
ingredients of the QCD factorization approach for describing B J / ipK* and then we pro- 
ceed to compute vertex and hard spectator interactions. The ambiguity of the experimental 
determination of spin amplitude phases is addressed in Sec. III. Numerical calculations and 
results are presented in Sec. IV. Discussions and conclusions are shown in Sec. V. 

II. B J/i)K* IN QCD FACTORIZATION 
A. Factorization formula 

The general B ^ J / tpK* amplitude consists of three independent Lorentz scalars: 

A[B{p) J/tp{ej/^,pj/^)K*{eK*,PK*)] oc e*fi^eK,{ag^u + hp^p^ + ice^uapp'}/^PK*)^ (2-1) 

where e''^^^ = +1 in our convention, the coefficient c corresponds to the P-wave amplitude, 
and a, h to the mixture of S- and Z^-wave amplitudes. Three helicity amplitudes can be 
constructed as [| 



*For B — > J l%\)K* decay the transverse amplitudes are given by li± = —a it uibPcC- 
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2mj/^mK* 
H± = a± rriBPc c, 



B 



m 



J/V 



(2.2) 



where Pc is the cm. momentum of the vector meson in the B rest frame. If the final-state two 
vector mesons are both hght as in charmless B ^ V1V2 decays with Vi being a recoiled meson 
and V2 an ejected one, it is expected that |-f^oP > l-^+P > l-f^-P owing to the argument that 
the amplitude if+ is suppressed by a factor of v^m2/mB as one of the quark helicities in 
V2 has to be flipped, while the H_ amplitude is subject to a further chirality suppression of 
order mi/niB 0. However, for B J /ijjK* decay, \/2mji^/mB is of order unity and hence 
in practice Hj^ and Hq can be comparable. 

Note that the polarized decay amplitudes can be expressed in several different but equiv- 
alent bases. For example, the helicity amplitudes can be related to the spin amplitudes in 
the transversity basis (Aq, Ay, A_|_) defined in terms of the linear polarization of the vector 
mesons, or to the partial-wave amplitudes (5, P, D) via: 



where we have followed the sign convention of [|T^]. The decay rate reads 



T{B J/i)K*) 



^K,v;:i^(|Ao|2 + |A^|2 + |Ai||2), 



Pc ,Gf, 



Pc ,Gf, 



B 



V2 



V,tV:,\'{\S\'+\P\'+\D\'). 



The effective Hamiltonian relevant for B —>■ J/ipK* has the form 



G, 
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Cl(/U)Oi(;U) + C2(yU)02(/i) - V^foV^^* ^ Ci(/i)Oi(/i) > + h.c. 



1=3 



where 



(2.3) 



(2.4) 



(2.5) 



Oi = icb)^._^{sc)^_^ 

o 



3(5) = ls6)^_^^(gV)v-A(v+A)) 
q' 



7(9) = 2^*^)^-aZ! V(5V)v+a(v-a)' 



O2 = {sb)^_^{cc)^_^, 

Oi(G) = {sab^)y_^^iq'f3qa)y_A{y+A)^ 

^8(10) = 



Sabfs)y_^J2(^q'Wc 

q' 



(2.6) 



aJv+A 



(v-a)' 
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with being the QCD penguin operators, O7-O10 the electroweak penguin operators, 

and (^ig'2)v±A = ^i7Ai(l i 75)^2 • Under factorization, the decay amphtude oi B ^ J/ifjK* 
reads 

A{B ^ J/i,K*) = ^V,,V:Xa2 + 03 + as + 07 + a,)X^''''' (2.7) 

where 



• (2.8) 



Note that for B J/ipK* decay, the factorizable amphtude X^^^ ''^^'^^ = 

{J/ip\{cc)y^^\0){K*\{sb)y_^\B) is the same as ( p. 81 ) except that the last term proportional 
to ie^uap has a positive sign. Comparing ( p. 81) with (|2.2|) leads to the helicity amplitudes 



Ho = -a{J/ipK*)ho, H± = d{J/ipK*)h±, (2.9) 

where a{J/'ijjK*) = 02+03+05+07+09. Note that the helicity amplitudes H± in B ^ J/ipK* 
are precisely the ones in 5 — > J/ipK* decays. Hence, in the factorization approach one 
has \H_\ > \H^\ for the former and > \H_\ for the latter. This is consistent with the 
picture that the s quark produced in the weak process b ccs in B J/ipK* has helicity 
— 1/2 in the zero quark mass limit. Therefore, the helicity of K* in B ^ J/ipK* cannot be 
+ 1 and the corresponding helicity amplitude if+ vanishes in the chiral limit |Tl|. 



B. QCD factorization 

Under naive factorization, the coefficients Oj are given by 02i = C2i + -^021-1, a2i-i = 
C2i-i + -^C2i- Hence, al^^J /ipK*) = a2{J/ipK) for h = 0,+, —. In the present paper, we will 
compute nonfactorizable corrections to a2{J /ipK*). The effective parameters of entering 
into the helicity amplitudes Hq and H± are not the same. 

The QCD-improved factorization approach advocated recently in [|1^ allows us to com- 
pute the nonfactorizable corrections in the heavy quark limit since only hard interactions 
between the (-BVi) system and V2 survive in the mh 00 limit. Naive factorization is recov- 
ered in the heavy quark limit and to the zeroth order of QCD corrections. In this approach, 
the light-cone distribution amplitudes (LCDAs) play an essential role. The LCDAs of the 
vector meson are given by [0,0 



5 



{ViP,eMxh^q'iO)\0) = fvmv / d^e 



e ■ X 



{V{P,eMxh,^,q'mO) = -my i fv - 



nia + ma' 



my 



ifyi^P^i'^v PuXfj,) 



e ■ X 



{P-x 



(F(P,.)|g(x)g'(0)|0) = zl (fy - § !^^^±^ 

z V my 



e*-x)m''y I d^e'^'^-'^hPiO, (2.10) 



where = 0, ^ is the hght-cone momentum fraction of the quark q in the vector meson, fy 
and fy are vector and tensor decay constants, respectively, but the latter is scale dependent. 



In Eq. (gJ^), (0 and $^(0 are twist-2 DAs, while and are twist-3 ones. 

Since 



e ■ X 
P ■ X 



P^ = e{: + 



e ■ X m 



P-x 2P-X 



^x^ 



[2.11] 



it is clear that to order 0{my /m\) the approximated relation P'^ = eji* holds for a light 
vector meson, where ej^ is the polarization vector of a longitudinally (transversely) 
polarized vector meson. Also, to a good approximation one has = Py/my for a light 
vector meson like K*, Hence, P ■ e± = and Eq. (|2.10|) can be simplified for K* as 



{K*{P,eMxh,smO) = fK*mK* I' d^e^^^^ kii^flO + <^^?f ^"^(0 



{K*{P,eMxh,^,s{0)\0) = ^mK^fK*e,.^perP''x^ 
{K*{P,eMx)a,.smO) = -tfU^^Pu - el^P,) [' d^e'^^'-'^T iO, 

•J 

(ir*(P,e)|g»s(0)|0) = -\fK*mK^ e^^'^'^'/i^HO, 

/ Jo " 



(2.12) 



where h'{^) = dh{^)/d^ and we have neglected light quark masses and applied the relation 



£ ' X 2 ^ ' X 



{P-X 



-^my 



p-x 



{P^P, - P,P^) + (£^|,P, - £,||P^), 



(2.13) 



which vanishes for a light vector meson. From Eq. (|2.12|) we see that the twist-3 DA h^^ of 
K* does not make a contribution. 

In the heavy quark limit, the B meson wave function is given by 



(0|6„(x)g^(0)|P(p))U^=.,=o=-^[(|^ + mB)75]/3, / rfpe-^^^+^-[$f(p)+^_<l>^(p)],„, (2.14) 



with n_ = (1, 0, 0, —1) and the normalization conditions 



/ f/p$f(p) = l, / rfp<|.f(p) = 0. (2.15) 

JO Jo 

Likewise, to the leading order in l/rric, the J/ijj wave function has a similar expression 

X fd^e-'^^^^~[^(/'^{0 + 4-^2^^iOU (2.16) 



Since the J/ip meson is heavy, the use of the light-cone wave function for J ftp is problem- 
atic. The effects of higher twist wave functions have to be included and may not converge 
fast enough. Because the charmed quark in J ftp carries a momentum fraction of order 
~ rric/mj/^, the distribution amplitudes of J/ip vanish in the end point region. In the fol- 
lowing study we adopt $|| as the DA of the non-local vector current of J/ip rather than g^^^ 
as the DA of the e± component since the latter does not vanish at the end point. Hence, we 
will treat the J/ ip wave function on the same footing as the B meson. Comparing Eq. ( p.l6| ) 
with Eq. ( p.lOp we see that at the leading order in l/rric one has 

$i^''(0 = ^'/""iO = ^'/""iO, fj/^ = fj/^- (2.17) 

The inclusion of vertex-type corrections and hard spectator interaction in QCD factor- 
ization leads to 
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= C2 + 


Cl 

Nc 


a, Cp ^. 
+ — — Cl F^, 
All N, ' 






= C3 + 


C4 

Nc 






4 


= C5 + 


Cq 
Nc 




-12) 


4 


= C7 + 


C8 

Nc 




-12) 


4 


= Cg 


ClO 


as Cp . 





(2.18) 



where Cp = {N^ — l)/{2Nc) and the superscript h denotes the polarization of the vector 
mesons: /i = for helicity state, and h = ± for helicity ± ones. In the naive dimensional 
regularization (NDR) scheme for 75, F^ in Eq. ( p.l8|) has the form 

F'^ = -12 In - 18 + + f'^r, (2.19) 
nib 

where the hard scattering function fj arises from vertex corrections [see Figs. l(a)-l(d)] and 
fjj from the hard spectator interactions with a hard gluon exchange between the emitted 
vector meson and the spectator quark of the B meson, as depicted in Figs. l(e)-l(f). 
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•s b 



(a) 



(6) 




s b 




(d) 




q q 




(/) 



FIG. 1. Vertex and spectator corrections to ^ J/ipK* . 



C. Vertex corrections 



The calculation of vertex corrections in Fig. 1 is very similar to that m B ^ J/ipK 
decay and the detail can be found in [Q. In terms of the two hard kernels // and gi given by 



and 



2zi , , Inf 

^ + (3 - 20 



+ 



+ 



^(1-0 
1 



1-e 



2zi 



i-zi 1-^(1-0 [(1-^(1-0]% 

2^2^ \ ln(l -z)-m 



In z^ 



+ \?>{l- z) + 2zi 



1-^(1-0; 1-^(1-0 j 

In^ AzrlnzC, ln(l — 2;) — ivr 

Azr- 



g,= /'df *[;'*(?) 

JO " 



-4^ 



lne + 



1 , 2(1 + ^-22:0 \ , . 

zc, In z^ — tiT 



1-^(1-0 j' 



ln(l -z) + 



(2.20) 



1 



(1-zO' 



z^ 



[l-z{l-OV {l-z){l-zO' 



[l-z{l-0]\ 
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where r = fj^^mc/{fj/^mj/^), z = nij^^/m^, the first scattering function fj induced from 
vertex corrections has the form 



f^ = fj + gj{l 



ft = fi, (2.22) 



where 



(2.23) 

In writing Eqs. ( p.2C ) and ( 2.21 ) we have distinguished the contributions from and 
^"^^ for reader's convenience, though later we will apply Eq. ( p.lT] ). Also notice that we 
have applied the relation |^[[] 

r^4^l^ = 2(^y = 2e. (2.24) 

Three remarks are in order, (i) As shown in [Q, the transverse DA (^'^/^ contributes 
not only to the transverse amplitudes H± but also to the longitudinal amplitude Hq, and 
vice versa for the longitudinal DA ■ This occurs because J /if) is heavy: the coefficient 
in front of $|| in Eq. (|2.11| ) consists of not only the longitudinal polarization but also the 
transverse one. (ii) It is easily seen that in the zero J/ip mass limit, 

/o ^ J\^^J/^^^) l^gi^lne-S^Trj , (2.25) 

in agreement with [|12| for B vrvr, as it should be. (iii) The expression of Aq/A^ in Eq. 
( p.22| ) can be further simplified by applying equations of motion. Neglecting the mass of 



^It is known from HQET that below the file scale, where mc is the pole mass of the charmed quark, 
the vector and tensor currents receive the same anomalous dimensions; that is, fj^^ and fj/^rric 
scale as the same power. Up to the scale, rescales with a factor [as{mb)/as{mc)]*^^^^\ rUc 
with [as{mb)/as{rnc)]^/^, and the ratio of /Jy^//j/^ becomes [as(jnc)/asimb)f/^'' 2mc{mb)/mj/^ = 
(1.1 — 1.2) X 2mcimh)/mj/^, where b = (llA'^c — 2nj)/3 and mc{mh) is the running charmed quark 
mass at the nib scale. However, the scale factor [agijUc] / (y.s{'mbyf'^^^^^ ~ (1.1 — 1.2) is relatively 
small and can be neglected for our purposes. 
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light quarks, applying the equation sj^k*(1 — 75)^ = and sandwiching it between the K* 
and B states leads to the result: 

and hence A^^* {m^i^) / A^^\rn?ji^) = 1. Consequently, fj=fj + gi{l — z). 

D. Hard spectator interactions 

For hard spectator interactions, we write 

/// = fii(2) + ///(3), (2.27) 

where the subscript (...) denotes the twist dimension of the LCDA. To the leading-twist 
order, we obtain 

X p-r] + {p-2^ + r])z + Aez ^g) 

p{p ~ V + vz)[{p - Oip - v) + ivp-v^- pO^Y 

This can be further simplified by noting that p ~ C(AQCD/"^b) — >^ in the mi, ^ 00 limit. 
Hence, 



11(2) — ^TT 7 / = / ? / 

^ ' Nc Ois[p) tiQ Jo p Jo 4 JO T] 



where the z terms in the numerator cancel after the integration over ^ via Eq. ( 2.24| ). 
Likewise, for transversely polarization states, we find 



X 



/' d-pd^dv^f{-p)^'/^'{o^T{v) -krr^y (2.30) 

JO p?7"^(l — z) 



Note that the hard gluon exchange in the spectator diagrams is not as hard as in the vertex 
diagrams. Since the virtual gluon's momentum squared there is k"^ = {—pps + VPk*Y ~ 
—pfim\ ~ —Ahmb, where Ah is the hadronic scale ~ 500 MeV, we will set as ~ as{\/ Ahmb) 
in the spectator diagrams. The corresponding Wilson coefficients in the spectator diagrams 
are also evaluated at the ph = VAhnib scale. As for twist-3 contributions to hard spectator 
interactions, we find 

///(3) = 0, (2.31) 
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and 



Since asymptotically $^*(?7) = 677(1 — rf), the logarithmic divergence of the 77 integral 
in Eq. (|2.29 ) implies that the spectator interaction is dominated by soft gluon exchanges 



between the spectator quark and the charmed or anti-charmed quark of J/ip. Hence, QCD 

ftm- 



factorization breaks down even at the twist-2 level for ftuo-) ■ Thus we will treat the divergent 



integral as an unknown "model" parameter and write 



rf!5. + (2.33) 



with ph being a complex number whose phase may be caused by soft rescattering . Note 



that linear divergences are cancelled owing to the relation ( p.24| ). Needless to say, how to 



treat the unknown parameter is a major theoretical uncertainty in the QCD factorization 
approach. 



E. Distribution amplitudes 



If we apply the asymptotic form for the vector meson's LCDAs [T^ 



^l(x) = = g);^'{x) = 6x(l - x) 



(«)/ 



1 + (2x - 1)^ 



(2.34) 



it is easy to check that /7/(3) = 0. Since the scale relevant to hard spectator interactions 
is of order ph = V^ft^fc ~ 1-5 GeV, it is important to take into account the evolution of 
LCDAs from p = 00 down to the lower scale. The leading-twist LCDA $m can be expanded 
in terms of Gegenbauer polynomials C^''^ ||13|| : 



^m{x,p) = 6x(l - x) ( 1 + ^ a*^(/i)C2'f (2x - 1) 



(2.35) 



n=l 



where the Gegenbauer moments a„ are multiplicatively renormalized. To = 2 we have 



^^{x,p) = 6x{l — x) 
$^(x,yu) = Qx{l — x) 



l + 3ak + -4(5^2-1) 



(2.36) 



where ^ = 2x — 1. For twist-3 DAs we follow [|l^] to use| 



g±\x,fi) = 6x{l-x) ^ , , I , i - 

+ 65+[3a;(l — x) + {1 — x) ln(l — x) + xlnx] 
+ 65_ [(1 — x) ln(l — x) — X In x] , 



16 



')}(5e^-l) 



(2.37) 



+ 
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.11 



112 



15 
64 



C3 (3^3'' - 



(3 - 30^2 + 35^^ 



+ ^6+[2 + Inx + ln(l - x)] + ^^-[2^ + ln(l - x) 



Inxl 



where the Gegenbauer moments and couplings 773, <jj\''^ , 5+ _ for K* at the scale /i^ = 1 GeV^ 
and ^ = h GeV^ can be found in |T5|. It turns out that the end-point behavior of g^'^ for 
K* is substantially modified and is very different from that of the asymptotic form (see Fig. 
3 of 111). 



III. EXPERIMENTS 

The angular analysis of — > Jj-ipK*^ and —>■ J/ipK*^ has been carried out by 
CDF §, CLEO § and most recently by the B factories BaBar [§ and Belle §. The 
three polarized amplitudes are measured in the transversity basis with results summarized 
in Table IV. Experimental results are conventionally expressed in terms of spin amplitudes 
^o,±,|| normalized to unity, \Aq\'^ + + = 1. Since the measurement of interference 
terms in the angular distribution is limited to Re(74||74Q), Im(y4_L74Q) and Im(74_L74||), there 
exists a phase ambiguity: 



0_L ^ ±7r - 
0± - 011 ^ ±7r - (0_L - 0||). 

Take the BaBar measurement as an example : 

= -0.17 ±0.17, (/)|| = 2.50 ± 0.22 , ^ \H+\<\H^l 

where the phases are measured in radians. The other allowed solution is 

(/),= -2.97 ±0.17, = -2.50 ±0.22, ^ \HA>\H\ 



(3.1) 



(3.2) 



(3.3) 



•f-Note that there is a slight difference for the expressions of in [^] and |14] 

12 



As pointed out in |]TT|, the solution (|3.2|) indicates that Ay has a sign opposite to that of A± 
and hence < \H^\, in contradiction to what expected from factorization. Therefore, we 
will compare solution (|3.3| ) with the factorization approach. Obviously there is a 3-0" effect 
that 0|| is different from vr and this agrees with the CDF measurement. However, such an 
effect is not observed by Belle and CLEO (see Table IV). In Table IV we will only list those 
amplitude phases from solution (|3.3|) . 
The measured branching ratios are 

(13.7 ± 0.9 ± 1.1) X 10-^ BaBar [5] 
B{B+ J/i)K*+) = { (12.9 ± 0.8 ± 1.2) x 10"^ Belle || (3.4) 

(14.1 ± 2.3 ± 2.4) X 10-^ CLEO § 



and 



(12.4 ± 0.5 ± 0.9) X 10-4 BaBar [5] 
B{B^ J/i)K*^) = I (12.5 ± 0.6 ± 0.8) x 10"^ Belle g] . (3.5) 

;i3.2 lb 1.7 ± 1.7) X 10-4 CLEO ~ 



IV. NUMERICAL RESULTS 



To proceed we use the next-to-leading Wilson coefficients in the NDR scheme 



ci = 1.082, C2 = -0.185, C3 = 0.014, a = -0.035, C5 = 0.009, ce = -0.041, 

cy/a = -0.002, cs/a = 0.054, cg/a = -1.292, cio/a = 0.263, Cg = -0.143, (4.1) 



at /i = mfe(mfe) = 4.40 GeV for A^^ = 225 MeV taken from Table XXII of IH with a being 



an electromagnetic fine-structure coupling constant. For the decay constants, we use 
/k* = 221 MeV, fj/^ = 405 MeV, /b = 190 MeV, 



(4.2) 



and we will assume fy = fv for the tensor decay constant. For LCDAs we use those in 
Sec. II. E, and the B meson wave function 



^fip) = NBp'il-pyexp 



1 fpniBV 
'2 \ uob ) 



(4.3) 



with ujb = 0.25 GeV and Nb being a normalization constant. 

In the following study, we will consider eight distinct form-factor models: the Bauer- 
Stech-Wirbel (BSWI) model piPl. the modified BSW model (referred to as the BSWII 



model) IT9I, the relativistic light-front (LF) quark model the Neubert-Stech (NS) model 
pT| , the QCD sum rule calculation by Yang |2^, the Ball-Braun (BB) model based on the 
light-cone sum rule analysis the Melikhov-Stech (MS) model based on the constituent 
quark picture 0] and the Isgur-Wise scaling laws based on the SU(2) heavy quark symmetry 
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TABLE II. Form factors Af^* , and V^^* at = and = f^'^ji^ in various form-factor 

models. 





BSWI 


BSWII 


LF 


NS 


Yang 


BB 


MS 


YYK 


Af^*(0) 


0.33 


0.33 


0.26 


0.30 


0.18 


0.34 


0.36 


0.49 




0.45 


0.45 


0.37 


0.39 


0.24 


0.43 


0.43 


0.49 




0.33 


0.33 


0.24 


0.30 


0.17 


0.28 


0.32 


0.30 




0.46 


0.63 


0.43 


0.48 


0.31 


0.45 


0.50 


0.42 




0.37 


0.37 


0.35 


0.30 


0.21 


0.46 


0.44 


0.39 


^^^*(-^./^) 


0.55 


0.82 


0.42 


0.51 


0.40 


0.86 


0.77 


0.87 



(YYK) so that the form factor Ai is mostly flat, A2 is a monopole-type form factor and V is 
a dipole-type one [^. The values of the form factors Af^*, A^'^* and V^^* at = and 
= Tn?ji^ in various form-factor models are shown in Table II. 

Among the eight form-factor models, only a few of them are consistent with the lattice 
calculations at large g^, constraint from B (pK* at low and the constraint from heavy 
quark symmetry for the form-factor q"^ dependence. The BSWI model assumes a monopole 
behavior (i.e. n = 1) for all the form factors. However, this is not consistent with heavy 
quark symmetry for heavy-to-heavy transition. The BSWII model takes the BSW model 
results for the form factors at zero momentum transfer but makes a different ansatz for 
their dependence, namely a dipole behavior (i.e. n = 2) is assumed for the form factors 
Fi, Aq, A2, V, motivated by heavy quark symmetry, and a monopole dependence for Fq,Ai. 
However, the equality of the form factors Af^* and A^^* at = is ruled out by recent 
measurements of B (pK* decays ^ . Lattice calculations of V^^* , Aq^* and Af^* at large 
p4[ in conjunction with reasonable extrapolation to g^ = m'^j/^ indicate that V^^* {m^^^) 



is of order 0.70-0.80. 

The parameters a^{J/tpK*) defined by 

d^iJ/ipK*) = 4 + 4 + a'^ + at^ + a'^ (4.4) 

are calculated using Eq. ( ^.18| ) and their results are shown in Table III. Since the penguin 
parameters ^3579 are small, in practice we have ~ Note that 5° and are inde- 
pendent of the parameter pn introduced in Eq. (|2.33| ); that is, they are infrared safe. Since 
h- is quite small due to the compensation between the Af^* and V^^* terms and ///(s) is 
inversely proportional to h_, a~ becomes more sensitive than a+ to the form-factor model 
chosen. 

From the experimental measurement of spin amplitudes, it is ready to extract the param- 
eters in various form-factor models. We use the averaged decay rate T{B — > J/ipK*) = 
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TABLE III. The calculated parameters d (J/tpK*) {h = 0,+,—) for B J/ijjK* decay in 
QCD factorization using various form-factor models for B — K* transition. The experimental 
results for d^{J /ijjK*) are obtained using the averaged branching ratio oi B ^ J /il)K* measured 
by BaBar, Belle and CLEO in conjunction with the central values of the BaBar measurement for 
the spin amplitudes |^o,±,||P- Only the central values of Og^pt ^ne shown here. 







50 


1 " 1 cxpt 




a+ 


1 lexpt 


d 


1 1 expt 


BSWI 


0.11 


- i0.06 


0.19 


0.16 


- z0.05 


0.18 


-0.01 + i0.05 


0.06 


BSWII 


0.15 


- i0.06 


0.25 


0.14 


- z0.05 


0.15 


-0.07 + Z0.05 


0.14 


LF 


0.14 


- i0.06 


0.25 


0.19 


- i0.05 


0.23 


-0.02 + ^0.05 


0.07 


NS 


0.14 


- i0.06 


0.25 


0.18 


- z0.05 


0.20 


-0.03 + i0.05 


0.08 


Yang 


0.23 


- i0.06 


0.43 


0.25 


- i0.05 


0.30 


-0.16 + i0.05 


0.20 


BB 


0.12 


- i0.06 


0.20 


0.14 


- i0.05 


0.16 


-0.15 + i0.05 


0.23 


MS 


0.13 


- i0.06 


0.22 


0.14 


- z0.05 


0.16 


-0.07 + i0.05 


0.14 


YYK 


0.09 


- i0.06 


0.16 


0.13 


- i0.05 


0.15 


-0.06 + i0.05 


0.12 



(5.34 ± 0.23) X 10-^6 GeV obtained from Eqs. and and the central values of the 

spin amplitudes measured by BaBar as an illustration: 

|io|2 = 0.597 ± 0.028 ± 0.024, ji^p = 0.160 ± 0.032 ± 0.014, 

|i|| P = 0.243 ± 0.034 ± 0.017. (4.5) 

Then aP can be determined from Tl{B J/'XpK*) = T{B — > J/ipK*) x \Aq\'^ and likewise 
for a^. The results are shown in Table III. It is evident the "experimental" values of 
are polarization dependent: \aP\ > \d~^\ > \d~\, whereas the present QCD factorization 
calculation yields la"*"! > \aP\ > \d~\. 

Normalized spin amplitudes and their phases in B ^ J / tpK* decays calculated in various 
form-factor models using QCD factorization are exhibited in Table IV, where the unknown 
parameter in Eq. ( |2.33| ) is taken to be real and unity. For comparison, we also carry out 
the analysis in the partial wave basis as the phases of S, P and D partial wave amplitudes 
are the ones directly related to the long-range final state interactions. We see from the 
Tables that the predicted l^oP? I-^P; ^'^d branching ratios are too small, whereas \A±\'^ = 
|Pp is too large. It is also clear that a non-trivial phase 0|j deviated from — vr is seen 
in some form- factor models, but it is still too small compared to the BaBar measurement. 
Nevertheless, a large phase 0|| as implied by BaBar can be achieved by adjusting the phase of 
the complex parameter pn, but admittedly it is rather arbitrary. In other words, the present 
QCD factorization calculation cannot say something definite for the phase (f)\\. The partial 
wave decompositions 5", P and D corresponding to the relative orbital angular momentum 
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L = 0,1,2 between J/ip and K* uniquely determine the spin angular momentum. Our 
results are difficult to account for the observation : |Dp : |Pp ^ 3.5 : 1 : 1 from recent 
Babar and Belle measurements. 

TABLE IV. Normalized spin amplitudes and their phases (in radians) in B ^ J/ipK* decays 
calculated in various form-factor models using QCD factorization. The branching ratios given in 
the Table are for — > J/ipK*^. For comparison, experimental results form CDF, CLEO, BaBar 
and Belle are also exhibited. 











0± 


0|| 


i3(10-3) 


BSWI 


0.43 


0.33 


0.24 


-3.05 


-2.89 


0.76 


BSWII 


0.38 


0.36 


0.26 


3.13 


-3.12 


0.73 


LF 


0.41 


0.34 


0.25 


-3.09 


-2.95 


0.69 


NS 


0.40 


0.34 


0.25 


-3.10 


-2.99 


0.70 


Yang 


0.38 


0.36 


0.25 


-3.12 


-3.11 


0.64 


BB 


0.41 


0.34 


0.25 


-3.04 


-3.05 


0.77 


MS 


0.40 


0.35 


0.25 


-3.08 


-3.05 


0.75 


YYK 


0.44 


0.32 


0.23 


-2.99 


-2.95 


0.84 


CLEO @ 
CDF i 
BaBar @ 
Belle @ 


0.52 ±0.08 
0.59 ±0.06 
0.60 ± 0.04 
0.60 ± 0.05 


0.16 ±0.09 


U-J^-J-O.ll 

0.16 ±0.03 
0.19 ±0.06 


0.32 ±0.12 
0.28 ±0.12 
0.24 ±0.04 
0.21 ±0.08 


-3.03 ± 0.46 
-2.58 ±0.54 
-2.97 ±0.17 
-3.15 ±0.21 


-3.00 ±0.37 
-2.20 ±0.47 
-2.50 ±0.22 
-2.86 ±0.25 


1.41 ±0.31 

1.37 ±0.14 
1.29 ±0.14 



There are several major theoretical uncertainties in the calculation: B — K* form factors, 
the twist-3 LCDAs of K* at the scale fih and the infrared divergences occurred in twist-2 
and twist-3 contributions. It has been advocated that Sudakov form factor suppression may 
alleviate the soft divergence [^. Hence, we have studied Sudakov effects explicitly and the 
detailed results will be presented in a future publication. When partons in the meson carry 
the transverse momentum through the exchange of gluons, the Sudakov suppression effect 
will be naturally generated due to large double logarithms exp which win 

suppress the long-distance contributions in the small k± region and give a sizable average 
(^1) ~ Am^, where A = ms — mb. This can resolve the singularity problem occurring 
at the end point. Basically, there is no Sudakov suppression in the vertex correction since 
the end-point singularity in the hard kernel is cancelled in the convolution. However, for 
the hard spectator interaction, we can have large Sudakov suppression effects at the end 
point since there are sizable {kjj contributions in the propagators. Especially, the end-point 
singularities without k± do not compensate in the twist-3 contributions. We find that 
is suppressed whereas is enhanced by the Sudakov effect and conclude that Sudakov 
suppression cannot help to solve the discrepancy between theory and experiment. 



TABLE V. Normalized partial wave amplitudes and their phases (in radians) in B ^ J /ifjK* 
decays calculated in various form-factor models using QCD factorization and fitted from the data, 
where (j)p = avg{PS*), (pD = arg(DS'*) and there exists a phase ambiguity: cpD — > —4'D and 
(f>P — > zbvr — (j)p. 





1*^ 1 


|P|2 


|Z)|2 




(7) 


BSWI 


0.60 


0.33 


0.07 


-0.04 


2.75 


BSWII 


0.60 


0.36 


0.04 


0.02 


3.10 


LF 


0.60 


0.34 


0.06 


-0.05 


2.80 


NS 


0.60 


0.34 


0.06 


-0.05 


2.86 


Yang 


0.59 


0.36 


0.05 


0.002 


3.07 


BB 


0.60 


0.34 


0.06 


0.05 


2.99 


MS 


0.60 


0.35 


0.05 


0.01 


2.97 


YYK 


0.60 


0.32 


0.08 


0.05 


2.85 


CLEO II 


0.77 ±0.19 


0.16 ±0.09 


0.07 ±0.03 


0.04 ± 0.59 


2.9 ±0.59 


CDF i 


0.61 ±0.34 




'-'•^•5-0.11 


0.26 ±0.20 


0.10 ±0.34 


2.17 ±0.34 


BaBar g 


0.65 ±0.13 


0.16 ±0.03 


0.19 ±0.10 


-0.13 ±0.21 


2.44 ±0.21 


Belle @ 


0.66 ±0.14 


0.19 ±0.06 


0.15 ±0.03 


-0.14 ±0.29 


2.80 ±0.29 



V. DISCUSSIONS AND CONCLUSIONS 

The hadronic decay B —* J/ipK* is analyzed within the framework of QCD factorization. 
The spin amplitudes Aq, Ay and Aj_ in the transversity basis and their relative phases are 
studied using various different form-factor models for B — K* transition. The effective 
parameters 02 for helicity /i = 0, ±, — states receive different nonfactorizable contributions 
and hence they are helicity dependent, contrary to naive factorization where a\ are universal 
and polarization independent. QCD factorization breaks down even at the twist-2 level 
for transverse hard spectator interactions. Although a nontrivial strong phase for the A\\ 
amplitude can be achieved by adjusting the phase of an infrared divergent contribution, the 
present QCD factorization calculation cannot say anything definite about the phase ^y. In 
QCD factorization we found that a\ and 02 are infrared safe. 

Unfortunately, our conclusion is somewhat negative, the longitudinal parameter a!^ cal- 
culated by QCD factorization which is of order 0.15 in magnitude is not large enough to 
account for the observed decay rates and the fraction of longitudinal polarization. In QCD 
factorization, the ratio R of vector meson to pseudoscalar production is close to unity with 
large uncertainties arising from the chirally enhanced and infrared sensitive contributions to 



B J/il)K (In the naive factorization approach, R ranges from 1.3 to 4.2 |g6|, but 
it is difficult to account for R, Tl/T, and |Pp simultaneously.) This is mainly ascribed to 
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the smallness of ag. It is instructive to compare a\{J /ipK*) in B ^ J/ipK* decay with 
a2{J/tpK) in B ^ J/tpK. It is found in § that a2{J/i)K) = 0.19l|j;}^ for \ph\ < 1 and 
that twist-2 as well as twist-3 hard spectator interactions are equally important. As for 
a2{J/i'K*), it is dominated by twist-2 hard spectator interactions. We have studied Su- 
dakov form factor suppression on end-point singularities and found that it does not help to 
solve the discrepancy between theory and experiment. 

Since the predicted in QCD factorization is too small compared to experiment, one 
may explore other effects that have not been studied. One possibility is that soft final-state 
interactions (FSIs) may enhance substantially [^]. A recent observation of 5^ D^*^*)-/?" 
decay by Belle p8[ and CLEO |29[ indicates a2{D^*^7r) ~ 0.40 - 0.55 much larger than 



the naive value of order 0.25. It is thus conceivable that some sort of inelastic FSIs could 
make substantial nonperturbative contributions to a^. The other possibility arises from the 
gluon component in the K* wave function. Consider the diagram in which one of outgoing 
charmed quarks emits a hard gluon before they form the J/ ip meson and the gluon fragments 
into a parton of the K* meson. Neglecting the charmed quark mass, because the charmed 
quark's helicity is conserved in the strong interaction, this gluon has zero helicity, i.e., it is 
longitudinally polarized. Following the same argument right after Eq. (|2.2|) , the hybrid K* 
will make a contribution to Hq and Although this amplitude is suppressed by order 
of Aqcd/'^6 owing to the presence of an additional propagator compared to the leading 
diagram, it is enhanced by the large Wilson coefficient Ci and hence cannot be ignored. A 
similar mechanism can also give a contribution to the B ^ J / ipK mode but it is difficult to 
make a quantitative estimate since the chirally enhanced twist-3 contribution is still quite 
uncertain. Good candidates to search for evidence of this effect are B , p^oj.ujuj. 

Without taking into account the hard gluon emission, the branching ratios of these decays 
which are color suppressed and dominated hy h ^ d penguin contributions are of order 
10"^ iO|,il],ll- Nevertheless, they can receive large contributions, proportional to ci at 
the amplitude level, from the hard gluon emission mechanism so that the branching ratios 
become 10~^ ~ 10~^. 



Note added: We learned the paper by X.S. Nguyen and X.Y. Pham (NP) ( |hep-ph/011028i , 
v2) in which a similar analysis in QCD factorization was carried out. However, their results 
differ from ours in some aspects: (i) There are some discrepancies between Eqs. ( |2.20[ - p?22D 
in the present paper and Eqs. (36) and (37) of NP. Also the expression of Ffj given by Eq. 
(39) of NP originally derived in |^ is valid only for two light vector mesons in the final state. 
It will get some modifications for heavy J/ip. It should be stressed that Eq. (28) adopted 
by NP for describing LCDAs works only for a light vector meson, but not for a heavy meson 
like J (ii) For hard spectator interactions, we have considered contributions from leading 
wave functions of B and J /if) and twist-3 DAs of K* [see Eq. ( |2.32D ], which are absent 



18 



in NP. Also we have taken into account the relevant scale /x/j = v^A/^m^ for hard spectator 
interactions, (iii) Unlike NP we did not consider the higher twist expansion for the J/ ip wave 
function. The twist expansion of LCDAs is applicable for light mesons but it is problematic 
for heavy mesons hke J/i/j. Note that although twist-3 J/t/j contributions to hard spectator 
interactions were considered by NP, they did not consistently compute the twist-3 effects of 
J/ijj in vertex corrections. 
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